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We compute the magnetic response of hollow semimetal cylinders and rings to the presence of an 
axial Aharonov-Bohm magnetic flux, in the absence of interactions. We predict nullification of the 
Aharonov-Bohm effect for a class of dispersion laws that includes "non-relativistic" dispersion and 
demonstrate that at zero flux the ground-state of a very short "armchair" graphene tube will exhibit 
a ferromagnetic broken symmetry. We also compute the diamagnetic response of bulk semimetals 
to the presence of a uniform magnetic field, specifically predicting that the susceptibility has a 
logarithmic dependence on the size of the sample. 



PACS numbers: 72.80.Vp, 73.23.-b, 73.23Ad, 73.63.Fg 

Ajiki and Ando (AA) 0] and Kane and Mele 2] have 
observed that the long-wavelength low-energy dynamics 
of the electrons of graphene Q when confined to a cylin- 
drical surface is described by the two-dimensional mass- 
less Dirac equation, and that the effects of the tube size 
and its chirality can be represented by a fictitious vec- 
tor potential. If additionally there is an axial Aharonov- 
Bohm (AB) flux $ present, then the energy eigenvalues 
of the Dirac equation are given by [l[ 
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where the overall upper and lower signs refer to the con- 
duction and valence bands (respectively), 7 is the Fermi 
velocity vp times fi, q z is the wave vector in the axial 
direction, W is the cylinder circumference, n — 0, ±1, ... 
is the azimuthal quantum number, and <f> = <&/ 4>o is the 
dimensionless AB flux measured in units of the flux quan- 
tum $0 = hc/e. The fictitious flux parameter a combines 
the effects of winding and curvature, and provides a clas- 
sification of nanotubes [1JJ2J . It has opposite signs in the 
K and K' Dirac valleys |l| as indicated by the lower and 
upper signs in front of a in Eq.([l]). This flux is not due 
to a physical magnetic field and for (f> = the system is 
time-reversal symmetric. All nanotubes can be classified 
as being semimetallic "armchair" (a — 0), insulating (a 
close to 1/3), or semiconducting (a -C 1) 

Since the integer part of a or <f> can be absorbed into 
the definition of the azimuthal quantum number n, any 
physical property is a periodic function of a or <f> with 
unit period. As an example of such a property, AA |l| cal- 
culated the AB magnetic response of an undoped cylin- 
der, which is a conceptually interesting problem because 
at zero temperature the free carriers are absent and the 
effect is entirely due to the electrons of the filled Dirac 
sea. Additionally, graphene, while being a semimetal, 
represents a marginal case between normal metals (where 
AB magnetism is experimentally interpretable in terms 
of persistent currents [4]) and insulators (where the effect 
is expected to be suppressed due to the band gap). For a 



non-chiral a = tube of length L, the valence electrons 
of given spin and belonging to the K- valley contribute to 
the ground-state energy the quantity 



£ 



(cyl) 

K 



(<t>) 



' E 

n=— oo 



Ldq z 
2tt 



2ir\ , 

w 1 



(2) 

The corresponding magnetic moment A4^- yl \(j>) and dif- 
ferential susceptibility XK^ifi) are determined by differ- 
entiation of Eq.([2]): 
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For a generic tube the ground-state energy (accounting 
for the spin degeneracy and including the contributions 
of both valleys) can be written in terms of Eg as 

£^ vl \<j),a) = 2{£ 
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and similar relationships hold for the total magnetic mo- 
ment and susceptibility, with £ replaced by A4 or x- Thus 
by computing only one of the if -functions as a function 
of the flux <j) we can understand the general problem. 

In what follows we will be also interested in the one- 
dimensional (ring) version of the same problem that is 
obtained by setting q z = in the spectrum Eq. (TT]) . Then 
instead of Eq. © the ground-state energy is given by 
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and relationships analogous to those of Eq.Q hold for 
the magnetic moment and susceptibility. 

At low energies the " relativistic" dispersion law ^ 
emerges in a variety of physical systems and so the 
problem of the AB magnetism of a semimetal goes be- 
yond graphene. We will give a comprehensive treatment 
of the phenomenon by employing the zeta function regu- 
larization method Q that finds wide applications in cal- 
culations of the Casimir effect. The flexibility and gen- 
erality of the technique will allow us not only to solve 
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the above-mentioned "relativistic" versions of the prob- 
lem but also, at no extra cost, discuss systems having 
more general dispersion laws. As a by-product we will 
also consider the Landau diamagnetism in semimetals. 

Direct inspection of Eqs.([2]) and ([5]) shows that they 
are divergent. The divergences are fictitious because the 
expression for the spectrum (fTJ) is only applicable at low 
energy; furthermore, the sum and integral should only 
be over wavevectors within the first Brillouin zone. AA 
[H treated this problem by introducing a cutoff function 
into Eq.@ which allowed them to carry out a numeri- 
cal calculation of the magnetic moment. From this they 
identified a cutoff-independent part which they argued 
captured the low energy part (Jf} of the true spectrum. 
A compact derivation of AA's result will be given below 
as a special case of a more general theory. 

Our calculation follows the analysis of a similar prob- 
lem . We begin by defining the spectral zeta functions 
for the cylinder 



A c v l ) 

SM 



n= — oo 

and the ring 



I.) ( " + ' 



M 



Jring) 



E 



(n + < 



(6) 



(7) 



versions of the problem. Here M confers a gap to the 
fermion spectrum, which will be put to zero later, and s 
is a parameter. For M finite and s positive and suf- 
ficiently large the expressions ^ and ([7]) are conver- 
gent and can be explicitly evaluated. The outcome will 
be analytically continued to the physically relevant sit- 
uation of M = and s = — 1. This procedure ex- 
tracts a cutoff-independent AB piece of the energies ([2]) 



and © via S^ yl, (4>) = -jL^ (-1) and £™(<£) 

"TCo"™ 9 '' ("1) which is all we need to compute magnetic 
properties ([3]). The contribution that is dropped rep- 
resents the cutoff-dependent ground-state energy of the 
filled Dirac sea in the absence of the AB flux; its value 
can be obtained by replacing the summation in Eqs.@ 
and ([5]) by an integration. 

The spectral zeta functions (JSJ) and ([7]) contain the so- 
lutions to other problems, too. Forty years ago Abrikosov 
and Beneslavskii [8[ demonstrated that (in three dimen- 
sions) crystal symmetry permits both linear s = — 1 (like 
in graphene) and parabolic ("non-relativistic") s = — 2 
touching of the valence and conduction bands. The latter 
parallels parabolic dispersion law found in unbiased bi- 
layer graphene Q if the interactions are neglected; more 
generally, s = — v describes a rhombohedral multilayer 
composed of v graphene monolayers [3, 0]- Abrikosov 
and Beneslavskii additionally investigated the role of 
Coulomb interactions whose effect, like in graphene [loj . 



was shown to be fairly weak in the case of a linear spec- 
trum. However, Coulomb interactions have a dramatic 
consequence for the case of a parabolic spectrum, s = — 2, 
where a breakdown of single-particle description was pre- 
dicted Q. The situation in bilayer graphene is simi- 



lar where recent experimental and theoretical work 11 1 
found that interactions can lead to a reconstruction of 
the ground state. Additionally, the s = —2 case war- 
rants special attention, because (i) the parabolic disper- 
sion plays an important role in an explanation of the un- 
conventional quantum Hall effect [9] and universal con- 
ductivity [12j in bilayer graphene, and (ii) it separates 
the regimes where the density of states is non-singular 
(for — s < 2) vs. singular (for — s > 2). 

The spectral zeta functions ^ and (UJ can be calcu- 
lated by using the identity 13] 

cospxdx f p \^r- K s^_i(pa) 

o (^W = ^) -rt|T' ns>0 

(8) 

where T(z) is the Gamma function and K^z) is the 
(MacDonald) modified Bessel function. For 5Rs > 1 this 
permits integration of ([6]) over q z , leading to 
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which relates the cylinder © and ring (J7J spectral zeta 
functions. Thus knowledge of the full s dependence of one 
of them describes both cases [13] • For example, the AB 
magnetism of a graphene ring (described by Cm=o'( — 1)) 
can be inferred via © from the zeta function for the 

cylinder (as Cm=o( ))- 

When JJs > 1, the spectral zeta function for the ring 
([Jj can be computed by employing the Poisson summa- 
tion formula. With the aid of (151) we find 
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Combined with (|9]) this provides us with an expression 
for the cylinder spectral zeta function 
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valid for 5Rs > 2. These expressions can be analytically 
continued into the 3?s < 1 and 5Rs < 2 regions, respec- 
tively, and the M — > limit can be taken which leads to 
our main results 
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As a first application of Eq. (|T2[) we consider the s = — 1 
case, which describes a cylinder with a linear dispersion 
law ([1]). The part of the ground-state energy that de- 
pends on the AB flux will be given by 
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The AB flux controls both the magnitude and sign of 
the result; the energy has maxima at <j) integer and min- 
ima at 4> half-odd integer. The magnetic moment and 
susceptibility of the cylinder follow from Eqs.© as: 
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where the last two representations, valid in the < (f> < 1 
range, should be periodically continued for all other <fi. 
We observe that both quantities are proportional to the 
cylinder length L; the magnetic moment is independent 
of the circumference W, while the susceptibility is loga- 
rithmically divergent at <fi integer. This translates into a 
weak non-analyticity at integer values of <p for the energy 
(fT4")l and magnetic moment (|T5|) . with the latter vanish- 
ing both at integer and half-odd integer cj). Eqs. (fl"5l) and 
(Jin]), combined with Eqs.@ for the total magnetic mo- 
ment and susceptibility, reproduce the AA results [l|. 

As an application of Eq.(fT3"|) we consider the s — — 1 
case which will describe a ring with linear dispersion law. 
The AB piece of the ground-state energy will be given by 
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where the last representation is valid for the range < 
4> < 1 and should be periodically continued for all other 
4>. As in the case of the cylinder, the magnitude and sign 



can be controlled by the AB flux <f>. The maxima of (fl7|) 
are located at 4> integer while the minima lie at half-odd 
integer </>. The magnetic moment then follows as 
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It is proportional to the circumference W, varies between 
jW/2&o and — jW/2$> , and has discontinuities at <fi in- 
teger. The susceptibility is constant and diamagnetic, 
except for positive delta-function peaks at </> integer. 

A remarkable feature of the ring geometry is that it 
displays a " ferromagnetic" broken symmetry at zero flux: 
the magnetic moment can be of either sign, depending on 
the history. This was already implied by Eq.([S]), which 
has discontinuous derivative at every integer value of 4>. 
A possible experimental realization of such a ring rep- 
resents an "armchair" graphene cylinder (a = case of 
Eqs.([T]) and ((H)) whose length is much smaller than its 
circumference, L <^ W . Such cylinders are not yet ex- 
perimentally available but hopefully peculiarity of their 
ground state would stimulate efforts to produce them. 

The case of the parabolic dispersion law s = — 2 holds 
a surprise: the cylinder and the ring spectral zeta func- 
tions (|T2"j) and (fT"3"|) (as well as the corresponding magnetic 
moments) vanish at s = —2. In fact, this remains true 
for any even v = — s because this is where the Gamma 
function in the denominators of Eqs. ([T2| and (fT3|) have 
poles. We thus conclude that the AB effect does not 
exist for a cylinder or ring with an even-layer rhombo- 
hedral graphene wall. Inspection of Eqs. flTUl) and (fTTj) 
shows that for even v = — s the AB effect is also identi- 
cally zero in the presence of a gap which covers the case of 
a dielectric or a filled band of a metal (s = —2, M 7^ 0). 

The AB magnetism is physically closely related to the 
Landau diamagnetism as the latter is also due to cur- 
rents circulating along the surface of the sample Q . The 
mathematical description of the two effects is also very 
similar. Indeed, for a linear dispersion law the energy 
eigenvalues of the three-dimensional Dirac equation are 
given by (compare with (P)) 
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where q z is parallel to the magnetic field H . Each of 
the Landau levels (fl"9)) labeled by n has a degeneracy 
eHA/2Trhc, where A is the cross-sectional area of the 
sample perpendicular to the direction of the field. In a 
semimetal, all negative energy Landau levels are filled 
while the n = state (shared between the valence and 
conduction bands) is half-filled. With this in mind the 
ground-state energy can be written as (compare with @) 
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where L is the height of the sample in the direction of 
the magnetic field. For a two-dimensional semimetal in a 
magnetic field perpendicular to the plane of the sample 
we would instead write (compare with ([5])) 
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As in the AB case, let us define the spectral zeta functions 
(compare with © and for M = 0) 
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Comparing Eqs.([7]) and (J53J) we notice that the latter 
can be analytically continued into the physically inter- 
esting region of negative s by setting = and replacing 
2n/W 2eH/hc and s s/2 in Eq.^: 
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where £(z) is the Riemann zeta function. The part of 
the energy (|2ip per unit area dependent on the magnetic 
field is then given by 



where, to give a better idea of the magnitude of the effect, 
we substituted 7 = Hvf ■ The magnetization A^*- 3d ' and 
magnetic susceptibility ^ 3d ) then follow as 
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The last equation erroneously predicts that as L — > 00, 
the susceptibility drops below the ideal diamagnetic limit 
of — 1/47T. This means that a more careful treatment 
is needed that distinguishes between the external mag- 
netic field H and the magnetic induction B representing 
the field experienced by the electrons of the substance 
18]. This can be accomplished by replacing H with B in 
Eq. (l2l?)) which no longer gives M(H); the latter depen- 
dence can be found from the equation H = B — AttM.{B). 
As a result correct version of Eqs. (|29|) would read 
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where we assumed that a cylindrical sample is placed in 
an external axial magnetic field H. We now see that in 
the thermodynamic limit L — » 00 the susceptibility x 
approaches — 1/47T, i.e. the bulk semimetal is an ideal 
diamagnet. In practice, however, we have |x| ^ 1 and 
Eqs. ([29| are adequate, as the amplitude of the logarithm 
in (|29|) is of the order 10 -6 ; astronomically large sample 
sizes would be required to observe |x| ~ 1/47T. 
This work was supported by the US AFOSR. 
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Multiplied by 4 (graphene's degeneracy factor) , this re- 
produces the result foreseen as early as 1989 fig . 

Since the spectral functions (|2"2"|) and (f2U)) satisfy the 
relationship (|9]), the latter combined with (p4|) provides 
us with the analytic continuation of (|22|) into the region 
of physically interesting s: 
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The physically relevant case now holds a surprise because 
at s = — 1 the spectral zeta function ([26]) has a pole: 
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This is a sign of a logarithmic cutoff dependence; the 
residue of the spectral function provides us with the am- 
plitude of the logarithm 17[ as indicated in the last step. 
Here b is of the order of the interparticle spacing. The 
magnetic piece of the energy (|2"0"|) per unit volume is then 
given by 
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